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ABSTRACT 

We investigate the non-Gaussian features of the IGM at redshift z 5 — 6 using Lya 
transmitted flux of quasar absorption spectra and cosmological hydrodynamic simulation of 
the concordance ACDM universe. We show that the neutral hydrogen mass density field and 
Lya transmitted flux fluctuations possess all the non-Gaussian features predicted by the log- 
Poisson hierarchy, which depends only on two dimensionless parameters (3 and 7, describing, 
respectively, the intermittence and singularity of the random fields. We find that the non- 
Gaussianity of the Lya transmitted flux of quasars from z ~ 4.9 to z = 6.3 can be well 
reconstructed by the hydrodynamical simulation samples. Although the Gunn-Peterson opti- 
cal depth and its variance underwent a significant evolution in the redshift range of 5 — 6, the 
intermittency measured by (3 is almost redshift-independent in this range. More interesting, 
the intermittency of quasar's absorption spectra on physical scales 0.1 — 1 h^^Mpc in redshift 
5 — 6 are found to be about the same as that on physical scales 1 — 10 h^^Mpc at redshifts 
2 — 4. Considering the Jeans length is less than 0.1 h^^Mpc at z ^ 5, and 1 h~^Mpc at 
z ~ 2, these results imply that the nonlinear evolution in high and low redshifts will lead 
the cosmic baryon fluid to a state similar to fully developed turbulence. The log-Poisson 
high order behavior of current high redshift data of quasar's spectrum can be explained by 
uniform UV background in the redshift range considered. We also studied the log-Poisson 
non-Gaussianity by considering inhomogeneous background. With several simplified models 
of inhomogeneous background, we found the effect of the inhomogeneous background on the 
log-Poisson non-Gaussianity is not larger than 1-sigma. 
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1 INTRODUCTION 

The cosmic density and velocity fields of the intergalactic 
medium(IGM) at low redshifts are highly non-Gaussian, while the 
temperature fluctuations of the cosmic microwave background ra- 
diation basically are Gaussian or at most weakly non-Gaussian. 
Thus, the non-Gaussian features of the IGM at low redshifts should 
mainly result from nonlinear dynamical processes. A well knowl- 
edge of these non-Gaussian features and their development history 
can fundamentally improve our understanding of the formation and 
evolution of the structures in the universe. For instance, the lognor- 
mal clustering model of the IGM (Bi, 1993; Bi & Devidson 1997; 
Liu et al 2006; Feng et al 2008) can explain well all the basic prop- 
erties of the Lya forests of quasar's absorption spectrum at various 
redshifts. However, high order statistics do reveal the deviation of 
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observation from the lognormal model (Jamkhedkar et al. 2003; Lu 
et al 2009). 

The aim of this paper is to study the non-Gaussian features 
of Lya transmission flux of quasar spectrum at high redshift. The 
non-Gaussianities of high redshift Lya transmission flux have been 
studied by various statistics, including the probability distribution 
function (PDF) of the flux (Fan et al. 2002; Becker et al. 2007), the 
distribution of the size of dark gaps (Songaila & Cowie 2002; Fan et 
al 2006), the largest peak width distribution (Gallerani et al 2007) 
and the distributions of the width of leaks (Liu et al 2007; Feng 
et al 2008). Most of these statistics are designed based on the fol- 
lowing observation facts: the high redshift quasar absorption spec- 
tra consist of complete absorption troughs separated by the spikes 
of transmitted flux. Obviously, these statistics are not suitable for 
comparing the non-Gaussian features at high redshift and low red- 
shift, as the absorption spectra of low redshift quasars consist of 
Lya forests and do not contain Gunn-Peterson troughs. Our focus 
is on the log-Poisson non-Gaussianity, which are characterized by 
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two dimensionless parameters that are available for both high and 
low redshift samples. Therefore, it can be a useful tool to study the 
redshift evolution of non-Gaussianity. 

More important, log-Poisson hierarchy is directly related to 
the dynamics of cosmic baryon fluid. With the cosmological hy- 
drodynamical simulation of the concordance ACDM model, it has 
been shown that, in the scale range from the onset of nonlinear 
evolution to the scale of dissipation, the velocity fields of cosmic 
baryon fluid at low redshift are extremely well described by the 
She-Leveque's scaling formula (He et al 2006), which was inferred 
from the log-Poisson hierarchical cascade (Dubrulle 1994; She & 
Wamire 1995; Benzi et al. 1996). The non-Gaussian behavior of 
the mass density field of the baryon fluid also can be well described 
by log-Poisson processes (Liu & Fang 2008). Recently, high order 
statistics of observed high resolution and high S/N Lya absorption 
spectra of quasars at redshift « ~ 2 - 3 are found to be well consis- 
tent with the non-Gaussian features predicted by the log-Poisson hi- 
erarchy (Lu et al 2009). It would be worth to study whether the sce- 
nario of self-similar log-Poisson hierarchical cascade is still hold 
at high redshifts. Whether the log-Poisson non-Gaussianity experi- 
ence a strong evolution in the redshift range 5-6? Can we explain 
the log-Poisson non-Gaussianity and its redshift dependence with 
the concordance ACDM model? 

A possible application of the non-Gaussianity of neutral hy- 
drogen distribution is to constrain the fluctuations of hydrogen- 
ionizing radiation background field. Considering the inhomogene- 
ity of UV background might be a source of the non-Gaussianity of 
ionized and neutral hydrogen, the non-Gaussianity of Lya transmit- 
ted flux would shed light on the evolution of UV background. The 
inhomogeneity of UV background at lower redshifts are negligible 
and a uniform ionizing background is an reasonable approximation. 
The situation becomes more complex and debatable when z > 5 
(Fan et al. 2002; Wyithe & Loeb 2005; Liu, et al. 2006; Fan at al. 
2006; Liu et al. 2007; Mesinger & Furlanetto, 2009; Furlanetto & 
Mesinger 2009). Therefore, it is also worth to study the effect of 
non-uniform UV backgroimd on the log-Poisson high order behav- 
ior at high redshifts. 

This paper is organized as follows. §2 gives the theoretical 
background. §3 presents the log-Poisson hierarchy of Lya trans- 
mitted flux fluctuations of observed samples of high redshift quasar 
absorption spectra. In §4, we describe the method used for produc- 
ing simulation samples, and show the log-Poisson non-Gaussianity 
of the neutral hydrogen component of cosmic baryon fluid. §5 com- 
pares the log-Poisson non-Gaussianities of observed data with sim- 
ulation samples of Lya transmitted flux. The effect of non-uniform 
UV background on non-Gaussianity is also studied in §5. Finally, 
conclusion and discussion are given in §6. 



2 THEORETICAL BACKGROUND 
2.1 Hierarchical clustering 

It has been recognized for a long time that the clustering process 
of cosmic matter probably is hierarchical (Peebles 1980). That is, 
the nonlinear dynamics of large scale structure formation can be 
characterized by the merging of holes from smaller to larger scales. 
An early model of hierarchical clustering assimies that the corre- 
lation functions of the density field satisfy a linked-pair relation, 
Cn = QnC"^^ (White, 1979), where is the nth irreducible cor- 
relation function with variable ^p(x) = p(x) — p, where p(x) and p 
are, respectively, the cosmic mass density field and its mean. How- 
ever, observation samples of transmitted flux of Lya forest do not 



support the linked-pair relation if the coefficients Q„ are assumed 
to be constants (e.g. Feng et al. 2001). 

The hierarchical merging has also been modeled by an addi- 
tive cascade rule (Cole and Kaiser, 1988). The basic step of the cas- 
cade rule is to assume that the a cell of mass M and spatial scale 
X will evolve into two cells, I and 2 with mass Mi = M + 5m 
and M2 = M — dm, on scale x/2, where 6m is a Gaussian ran- 
dom variable. However, central limit theorem shows that the field 
produced by an additive cascade process should be Gaussian. It 
cannot explain the non-Gaussian features of Lya transmitted flux 
(Pando et al. 1998). Moreover, this cascade needs a "very small ini- 
tial units" (Peacock, 1999) as the first generation halo of the hierar- 
chy, of which the mass is non-zero. This "initial units" is not com- 
patible with hydrodynamic equations with continuous variables. 

A proper model of the the hierarchical clustering should be 
randomly multiplicative, and infinitely divisible. For a randomly 
multiplicative cascade, the mass m„ in a cell at step n is related 
to step n — 1 by ?n„ = (1 ± where 5 is a random vari- 

able. Non-Gaussian features can be formed through randomly mul- 
tiplicative cascade processes, even if the original field is Gaussian. 
The infinite divisibility means that there is no finite "initial units" in 
the hierarchical merging. The "initial units" of the merging process 
can be infinitesimal. 

The log-Poisson hierarchy has all these desired properties. 
More important, the log-Poisson hierarchy actually was inferred 
from the invariance and symmetry of the Navier-Stokes equations 
and works well in explaining the high order behavior of fully devel- 
oped turbulence (Dubrulle, 1994; She & Waymire 1995; Leveque 
& She 1997). 



2.2 Log-Poisson hierarchy 

To measure the non-Gaussianity caused by log-Poisson hierarchy, it 
is better to use the variable (5pr ~ /9(x + r) — p(x),7- = |r|,butnot 
p(x -|- r) — p. For a statistically isotropic and homogeneous random 
field, one can just consider |5pT-|, as the distribution of positive and 
negative dpr is statistically symmetric. The basic statistical quan- 
tity is the structure function defined by 

Sp{r) = ilSprin, (1) 

where p is the order of statistics, and the average (...) is taken over 
the ensemble of density fields. The second-order structure function 
S2 = {|(5pr|^) as a function of r (scale) is actually the power spec- 
trum of the mass density field (Fang & Feng 2000). 

In scale-free range, the structure fimction should be a function 
of power law of r as 

Sp{r)ocr^^\ (2) 

where ^(p) is referred to intermittent exponent. For a Gaussian 
field, ^(p) is a linear function of p, but for a intermittent field, func- 
tion ^(p) is nonlinear. 

The log-Poisson hierarchy assumes that, in the scale-free 
range, the variables Spr on different scales r are related to each 
other by a statistically hierarchical relation as (Dubrulle 1994; She 
& Waymire 1995) 

5pr2 — W^rir2<5pri. (3) 

The factor Wrirs is a function of the ratio ri /r2 given by 
Wr,r, = rin/r^r, (4) 
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which describes how the fluctuation Spn on the larger scale ri 
relates to fluctuations 5pr2 on the smaller scale r2. In eq.(4), m is 
a Poisson random variable with the PDF 



P{m) = ex.p{-\r^r2)Kivi/r. 



(5) 



The random variable m can be considered as the steps of the evolu- 
tion from 5pr-^ to 5pr-2- To ensure the normaUzation {Wr^r^) = 1. 
where the average (...) is over m, the mean Anrg of the Poisson 
distribution is then 



7[Mri/r2)]/(l-/3). 



(6) 



The log-Poisson hierarchy contains both the spatial size and ampli- 
tude of the density fluctuations. This point is different from other 
hierarchical models, which consider only the size of hierarchical 
units. 

The log-Poisson hierarchy given by eq.(3) depends only on 
the ratio ri/r2, which is obviously scale invariant. The hierarchy is 
determined by two dimensionless positive parameters: /3 and 7, de- 
scribing, respectively, the intermittence and singularity of the ran- 
dom fields. Equation (3) relates 5pr on different scales by multi- 
plying a random factor W , which generally yields a non-Gaussian 
field even if the field originally is Gaussian (Pando et al. 1998). 

The cascade from scale ri to r2, and then to ra is identical to 
the cascade from n to rs. It is because = W^r2r3 = 

/3^(ri/r3)^, where A'^ is again a Poisson random variable with 
\r1r3 =Arir2+Ar2r3 = 7 [ln(r 1 /r3 )] / ( 1 — /3) . The log-Poissou 
hierarchy removes an arbitrariness in defining the steps of cascade 
from n to r2 or r2 to 7-3. Therefore, the log-Poisson hierarchy, sug- 
gested by eq.(3), is discrete in terms of the discrete random number 
m. However, the scale r is infinitely divisible. Namely, there is no 
lower limit on the difference ri — 7-2 . It can be infinitesimal, and the 
hierarchical process is of infinite divisibility. With the log-Poisson 
model eqs.(3)-(6), one can show that the intermittent exponent ^ (p) 
is given by (Liu & Fang 2008) 



?(p) = -7b-(l-/3")/(l-/3)]- 



(7) 



This is actually the SL scaling formula (She & Leveque 1994). 

From eq.(7), the power spectrum S2{r) = const is flat, or 
called white. However, the power spectrum of the initial Gaussian 
field of the cosmic matter generally is not white, but colored with 
power-law S2{t) cc r^^°. In this case, we should adjust the log- 
Poisson hierarchy eq.(3) by replacing Spn and 5pr2 with rfSpn 
and r2 5pr2 respectively. The intermittent exponent f (p) is 



(8) 



When parameter a = 0, eq.(8) will simplify to eq.(7). The non- 
Gaussian features of the field described by the log-Poisson hierar- 
chical clustering have been given in Liu & Fang (2008) and Lu et 
al (2009). 



3 LOG-POISSON NON-GAUSSIANITY OF OBSERVED 
SAMPLES 

3.1 Observed data 

Observational data used here consists of the spectra of 19 QSOs 
with redshifts from z =5.74 to 6.42 that compiled in Fan et al. 
(2006). The data have resolution of i? ~ 3000 - 4000 and are 
re-binned to a resolution R = 2600. The observed flux, /obs. is 
normalized with a power-law continuum /con oc The noise 

and continuum imcertainty of transmitted flux F = /obs //con is in 



the level of 0.018 ± 0.012 in the range z < 5.7, and 0.014 ± 0.008 
in the range z > 5.7 (Liu et al. 2007). For more details, we refer to 
Fan et al. (2006). Wavelength of the data covers roughly the rest- 
frame wavelength from 900 to 1350A. To avoid the mixing of Ly/3 
absorption and the proximity effect, only pixels have the rest frame 
wavelength > 1040 A and below the maximum Lya are used. 

In our analysis below, we use only the Lya transmitted flux 
in the redshift range from 4.7 to 6.3. There actually are 12 quasars 
available in the range z > 5.9 and reduce further to 4 in 6.1 < 
z < 6.3. We divide the redshift range from 4.9 to 6.3 into 7 bins by 
4.9 + n X 0.20 < « < 4.9 + (n + 1) x 0.20 where n = 0, 1, ... 6. 
More specifically, the redshift size of each bin is Az = 0.20. As 
a comparison, we also use sample in the redshift range from 4.7 to 
6.2, and divide it into 10 bins by 4.7 -|- n x 0.15 < z < 4.7 + 
(n + 1) X 0.15 where n = 0, 1, ... 9; each bin has a size of Az = 
0.15. All the transmission flux pixels in a given redshift bin form 
an ensemble. The numbers of pixels in different ensemble are not 
uniform. 

In order to compare with observation at moderate redshift 
we also analyze the high resolution and high signal to noise ra- 
tio Lya absorption spectra of 28 Keck High Resolution (HIRES) 
QSO (Kirkman & Tytler 1997). It is flie same as that used in Lu 
et al (2009). The details of the data set and its reduction have been 
described in Jamkhedkar et al. (2002, 2003). 

With these samples we calculate the optical depth t{z) = 
— In F(2;), and the fluctuation of the optical depth Stt = T{x+r) — 
t{x), where the spatial coordinates x and r are in physical scale. 
Since the variable Stt is given by difference between r(3; + r) and 
t{x), Stt is independent of fluctuations of t(x) on scales larger 
than r. Therefore, the variable Stt actually is insensitive to the con- 
tinuum used in the data reduction (Jamkhedkar et al. 2001). 

As has been shown in Lu et al (2009), the variable 5Tr is ap- 
proximately a measurement of the density fluctuation p,.- We can 
then study the log-Poisson non-Gaussianity with variable Stt. For 
instance, the structure function with the variable Stt is defined now 
by 



(9) 



To treat the unwanted data, including low S/N and bad pixels, 
we use the algorithm of wavelet denoising by threshold (Donoho 
1995, Jamkhedkar, et. al. 2003). This method is effective for pix- 
elated data. First, we calculate the wavelet scaling function coeffi- 
cients (SFCs) of both transmission flux field F(a;) and noise field 
n{x) with 

Eji = j ¥{x)4)j^i{x)dx, efi = J n{x)cl)j^i{x)dx. (10) 

where 4>j,i{x) is the scaling function of wavelet on scale j and 
at position 1. We then identify imwanted mode [j, I) by using the 
threshold condition 



I F / iVi , r 



(11) 



This condition flags all modes with S/N less than /. We skip all 
the flagged modes when doing statistics. To reduce the boundary 
effect of unwanted chunks, we also flag two models aroimd an un- 
wanted model. With this method, no rejoining and smoothing of 
the data are needed. The threshold / is given by the same way as 
Jamkhedkar et al (2003) and Lu et al (2009). 

A typical statistical quantity of log-Poisson non-Gaussianity 



Fp(r) = 5p+i(r)/5p(r). 



(12) 
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Figure 1. F4{r) given by observed sample of the Lya transmitted flux at redshift ranges 4.9 < z < 5.1 (left), and 5.7 < z < 5.9 (right). 



We test the effect of noise on F4{r) = S^i'r) / S4(r). We cal- 
culate F4{r) for data sets given by different threshold /. The /- 
dependence of F4{r) is shown in Figure 1. For the data set at 
4.9 < z < 5.1, the values of F4{r) are /-independent when 
log / > 0.3. In other words, the statistical results are stable with 
respect to threshold / ~ 2. We will use only data with S/N larger 
than 2. This threshold is larger than the error level given by Fan et 
al (2006). For sample set of 5.7 < z < 5.9 (right panel of Fig- 
ure 1), ^4(2;) are very weakly dependent on log/. It is because, 
in the redshift range z > 5.7, the high order statistics ^4(0;) are 
dominated by modes with high Stt, and therefore, it is insensitive 
to dropping modes with low Srr. However, the fewer the modes, 
the larger the variance of Poisson process. Therefore, we should 
consider the variance of Poisson process in our statistics. 

3.2 Redshift dependence of parameter /3 

We first study the /3-hierarchy predicted by the log-Poisson hierar- 
chy. It reads (Liu & Fang, 2008) 

In Fp+i{r)/F3{r)= 13 ]nFp{r)/F2{r). (13) 

Eq.(13) requires that for all r and p, ln[Fp+i(r)/F3(r)] vs. 
\n[Fp{r) / F2{r)] should be on a straight Une with slope /3. which 
is called /3-hierarchy. Eq.(13) does not contain parameters 7 and a. 
Figure 2 presents the /^-hierarchy of observed transmitted flux in 
4 redshift ranges z = 5.0, 5.4, 5.8 and 6.2. The statistical quan- 
tity Fp+i(r) are given by all available observed data, of which the 
physical scale r covers the range from ~ 0.1 to 1.5 h~^ Mpc, and 
the order parameter p increases from 1 to 2.5. All the distributions 
of Fp+i{r)/F3{r) vs. Fp{r) / F2{r) in Figure 2 can be well fitted 
with a straight Une. It shows that the Lya transmitted flux of real 
observation at high redshifts still satisfy the /3-hierarchy and indi- 
cates that the log-Poisson non-Gaussian features are significant in 
the considered redshift range. 

From the /3-hierarchy straight line, we calculate the parame- 
ter /3 in redshift ranges from z = 5.0 to 6.2. Figure 3(a) presents 
the mean of /3 in each redshift bin. The error bars are given by 
the variance of Poisson process, which generally are large than the 
variance of /3 in the given redshift bin. Figure 3(a) shows that the 
redshift evolution of /3 is rather weak in the range from 2 = 5 to 6. 
We also calculate the redshift-dependence of /3 with the same data. 



but the size of redshift bin is taken to be Az = 0.15. The result is 
also plotted in Figure 3(a), which shows the same redshift depen- 
dence as that of Az = 0.20. Therefore, the weak redshift-evolution 
of parameter B does not affected by Az. From now on, we will only 
give results with Az = 0.20. 

Figure 3(b) plots the redshift dependence of /3 for the Keck 
HIRES quasars spectra sample. The error bars are given by the 
maximum and minimum in each redshift interval. Comparing Fig- 
ure 3(a) and (b), we can conclude that the non-Gaussian parameter 
P evolves weakly from redshift 2 to 6. It should be pointed out that 
the physical scale r in Figure 3(a) (high redshift sample) is actually 
smaller than that of Figure 3(b) (low redshift sample). Figure 3 re- 
veals that the /3 non-Gaussianity on a small scale at high redshift 
are about equal to that on a large scale at low redshift. This result 
indicates that the turbulence state on scales 0.1-1 h^^Mpc at red- 
shift 2 ~ 5 would be the same as that on scales 1-10 h~^Mpc at 
redshift ~ 2. This is consist with the fact that the Jeans length at 
2 ~ 5 is less than 0.1 h^^Mpc and increases to ~ 1 h^^Mpc at 

2-2. 

Quasar's Lya absorption spectra at high redshift 2 ~ 5 — 6 
are significantly different from that at low redshift 2 ~ 2 — 3. The 
latter show Lya forests, while the former consist of complete ab- 
sorption troughs (Gunn-Peterson troughs) separated by tiny trans- 
parent regions. In other words, Lya transmitted tlux experience a 
strong evolution with redshift rising from 2~2 — 3to5 — 6. The 
low-order statistics of the Lya transmitted flux, such as the mean 
optical depth and its variance, also show strong redshift evolution. 
Therefore, the weak redshift-evolution of /3 is very interesting. It 
implies that the /3 non-Gaussian feature is mainly dependent on the 
nonlinear state of the fluid, but weakly dependent on the optical 
depth and its variance. 



3.3 Non-Gaussianity related to parameter 7 

We now turn to the 7-reIated non-Gaussianity, which is given by 
(Liu & Fang 2008) 

^r^jj^ = Kplnr + const, (14) 
and 
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Figure 2. The /3-hierarchy of observed sample of the Lya transmitted flux at redshift z = 5.0 (top-left) , 5.4 (top-right), 5.8 (bottom-left) and 6.2 (bottom- 
right). The physical scale r is in the range ~ 0.1 - 1.5 h"^ Mpc, and order p is from 1 to 2.5 . The error bars are given by the maximum and minimum of 
bootstrap re-samphng. 
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Figure 3. Redshift dependence of parameter ^ for samples of a.) Fan et al (2006) (left), and b.) the Keck HIREs quasars, Jamldiedkar et al.(2003) (right). 



X. = -7^^^^%f^. (15) 

It requires that for a given p, the relation of ln{{dTr)^^) / {{STr)^y 
and In r has to be a straight line with slope Kp. Since parameter /3 
is already determined by /3-hierarchy in last section, we can then 
figure out the parameter 7 from Kp. 

The result is presented in Figure 4. It shows the relation of 
lii{{5T^'') / ({6Tr)^y) vs. Inr for observational data at redshift 
ranges z = 5.0, 5.4, 5.8 and 6.2. The order parameter p is set to 2 
and 3, i.e. the statistics of eqs.(14) and (15) are of the order of 4 and 
6. The error bars are given by the maximum and minimum of each 
r. The In r-dependency of ln(((5T^^)/(((5Tr)^)'') approximately at 
each range can be given by a straight line. 



The slopes Kp of the straight lines of Figure 4 are listed in Ta- 
ble 1, and parameter 7 given by eq.(15) is also listed. In log-Poisson 
hierarchy, parameter 7 has to be independent of p and therefore, the 
values of 7 determined by Kp with different straight lines should 
be the same. Table 1 indeed confirms this point. We see that for a 
given redshift, the statistics of p = 2 and 3 yield the same 7 within 
their errors. Therefore, the high redshift Lya transmitted flux well 
fulfills the 7-related non-Gaussianity. 

A basic feature of log-Poisson model is to yield nonlinear 
terms of p, i.e. the term /3'' in eq.(7) and fi'^^ in eq.(15) (Frisch 
1995). The terms with linear p can also be given by other models. 
Since /3 < 1, the tests with P" and n > 6 don't give new test on 
the log-Poisson model. Therefore, the statistical order n generally 
is taken to be less than 6. 
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Figure 4. ln[(((5Tr)^^) /{(<5Tr)^)^] vs. r of observed sample of the Lya transmitted flux for redshifts at 5.0 (top-left), 5.4(top-right), 5.8(bottom-left), and 
6.2(bottom-right). p is taken to be 2 (bottom line) and 3 (top line). The solid lines are given by the least square fitting. The error bars are given by the maximum 
and minimum of bootstrap re-sampUng. 



Table 1. Parameter 7 at redshift z = 5.0 to 6.2 



5.0 



5.4 



5.8 



6.2 



Kp 

-0.46± 0.20 
-l.lOi 0.40 
-0.25±0.10 
-0.52± 0.30 
-0.24± 0.07 
-0.52±0.15 
-0.28±0.16 
-0.78± 0.35 



0.41 
0.43 
0.22: 
0.20 
0.21 
0.21 
0.27 
0.31 



+0.17 
-0.18 
-0.16 
-0.15 
+0.08 
0.07 
+0.12 
0.11 
+0.06 
-0.06 
-0.06 
-0.06 
+0.15 
0.15 
+0.13 
0.12 



4 HYDRODYNAMIC SIMULATION SAMPLES 

Although the mean optical depth and its variance of Lya transmit- 
ted flux of quasar's absorption spectrum underwent a strong evolu- 
tion at high redshift, the log-Poisson non-Gaussianity, as last sec- 
tion reveals, shows only weak dependence on redshift. Thus, an 
important question is whether the two aspects of the redshift evo- 
lutions can be conciliated within the concordance ACDM model. 
We study this problem with cosmological hydrodynamic simula- 
tion samples. 



4.1 Simulation samples of neutral hydrogen 

We first produce the samples of mass density, temperature and ve- 
locity fields of cosmic hydrogen with hydrodynamic simulation 
with the same code of Liu et al. (2007, 2008), which is based 
on Eulerian method for hydrodynamics with 5th order Weighted 
Essentially Non-oscillatory (WENO) finite difference scheme and 
particle mesh(PM) method for dark matter particles (Feng et al. 
2004). We use cosmological parameters given by the latest result of 
WMAP (Komatsu et al. 2009). We run simulations in a period box 



of side 25h~^ Mpc with 512^ grids and dark matter particles which 
have a mass resolution of 1.04 x IffM^. Assuming ionization 
equilibrium, atomic processes, including radiative cooling, heating 
and fraction of species, are modeled using the primordial compo- 
sition {X — 0.76, Y = 0.24) and formalism in the Appendix of 
Theuns et al. (1998), under optically thin approximation. Photoion- 
ization and photoheating are switched on after the UV backgrotmd 
is added at a = 11.0. 

Started atz = 99, a sample output at 2; = 11 first. Simulations 
with different UV histories are then performed from this snapshot 
at z = 11 and produce snapshots at redshifts from z = 6.5 to 
z = 4.9 at a interval of Az = 0.1. As we focus on the evolution 
at high redshifts, simulations are stopped at 2; = 4. From snapshot 
dumps, we produce the mass density field of hydrogen p(x), tem- 
perature field r(x), and velocity field v(x). Star formation is not 
included. However, the contribution of stars to the UV background 
is considered by fitting the redshift evolution of the UV background 
with mean optical depth and its variance of Lya transmitted flux 
field. This code has recently been used to produce samples to show 
the turbulence behavior at redshift as high as 2 ~ 4 (Zhu et al. 
2010). These samples would also be suitable to test the log-Poisson 
non-Gaussianity of the IGM at high redshifts. 

Although the Gunn-Peterson optical depth shows dramatic de- 
crease with redshift and abnormally large scatter at « ~ 6, it can 
still be fitted by models of a uniform ionizing background (Lidz 
et al. 2006; Liu et al. 2006, 2007; Mesinger & Furlanetto 2009). 
As current observation does not give a well knowledge of ioniz- 
ing source at high redshift, instead of the Haardt & Madau (2001) 
model of UV background history, we use a more general uniform 
UV background with the ionizing photons have a power-law spec- 
trum with index —1.0 and an normalized coefficient J21. The hy- 
drogen photoionization rate r_i2 is calculated using the fitting for- 
mula in Theuns et al. (1998) and then is used to calculating the 
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heating and cooling in the hydrodynamic simulation. The photoion- 
ization rate can be given by 

r_i2 =3.15J2i. (16) 

Based on the above assumption, the history of the UV back- 
ground in our simulation is given by the evolution of J21. There is 
no direct observation on this parameter. Considering the two typi- 
cal reionization history in references, one is the extended scenario 
and the other is phase transition, we use the following two redshift- 
dependent models of J21, 

J2i(^) = 5.0 X exp(-0.4752;). (17) 

J2i{z) = exp(-0.21^;^ -|- 1.5z - 3.0) -h 0.02exp(-7.0/z). (18) 

The J21 (z) of eq.(17) can be considered as a model of the extended 
reionization scenario. Eq(18) gives approximately the same evolu- 
tion as eq.(17) at z = 4.0 — 5.0, but drops about an order of mag- 
nitude from z = 5.0 to 2 = 7.0 and stays at very low level when 
a > 7. It is to mimic the reionization as a phase transition over 
redshift range from z = 5.0 to 2; = 7.0. The parameters used in 
eqs.(17) and (18) actually are determined by the fitting of simulated 
samples with observed mean optical depth and its variance of Lya 
transmitted flux field. Several simulations have been performed be- 
fore these parameters are selected. It is interesting to see that with 
the parameters of oqs.(17) and (18), the intensities of r_i2 given 
by eq.(16) in the redshift range 2; = 4 — 9 are just in between of 
the values given by 1.) Haardt & Madau, (2001) and 2.) proximity 
effect of quasars (Dall'Aglio et al. 2008, Gilmore et al. 2009 and 
references therein ). 

The last but not least, although the intensities of the UV 
background given by models (17) and (18) are very different at 
z = 5.0 to z = 7.0, the non-Gaussian features of the transmit- 
ted field of Lya are less affected. This is because the basic variable 
Stt = t{x + r) — t{x) is not sensitive to the change of the uniform 
ionizing background given by the models (17) and (18), especially 
when r is small. 

4.2 Intermittence of neutral hydrogen density field 

Before simulating the Lya transmitted flux, we analyze the log- 
Poisson behavior of the simulation samples of neutral hydrogen 
density field, nHi(x). We first calculate the /3-hierarchy of neutral 
hydrogen field. The results are presented in Figure 5, which con- 
tains all data points of fp+i (r) /^s (r ) and Fp (r) /F2 (r) with p = 
1, 1.5, 2, and 2.5, and available scale r. To estimate errors, we di- 
vide our 160 one-dimensional samples into 80 subsamples, each of 
which has 2 lines. The error bars are given by the scattering ranges 
of the 80 subsamples. Figure 5 shows that /3-hierarchy is also well 
hold for neutral hydrogen density field. The numbers of /3 are also 
shown in Figure 5. It is in the range 3.0 — 3.5. 

In log-Poisson hierarchy, the r-dependence of Fp{r) is given 
by (Liu et al 2008; Lu et al. 2009) 

Fp{r) ^ Ar-°'-''^''-^''\ (19) 

That is, the relation between In Fp{r) and In r should be a straight 
line with the slope -a + 7(1 — For Gaussian field we have 13 — 
1, and therefore, the slope is independent of p. For log-Poisson field 
(/? < 1), when p is large, the slope converges to — (a — 7). On the 
other hand, parameter a can be determined by the power spectrum 
of nmix) (§2.2). Thus, The \nFp{r)-r relation with larger p can 
be used to determine parameter 7. 

Figure 6 presents lnFj,(r) vs. Inr of the simulation samples 



of neutral hydrogen mass density field in the physical length scale 
range of 0.1 < r < 1.5 h~^ Mpc and orders of p = 0.5 x n with 
n = 1, 2. ..8. For all redshifts andp < 4, In Fj,(r) — Inr can indeed 
be approximately fitted by straight lines. These straight lines have 
different slopes. The higher the p, the more steeper the straight lines 
are. It implies that /3 < 1 and the density fields are non-Gaussian. 
When p > 3, the fitted straight fines are almost independent of p. 
It shows the convergent of the slope to -(a — 7) at high p. The 
parameters 7 found with slope of \aFp{r)-r lines are 0.70, 0.65, 
0.60 and 0.60. corresponding to redshifts z = 5.0, 5.4, 5.8 and 6.2 
respectively. 

From the parameters /3 and 7 given above, we plot the inter- 
mittent exponent ^(p) [eq.(8)] in Figure 7, in which data points are 
given by fitting structure functions of simulation samples to eq.(2) 
with p = 0.5 + 0.5n, n = 0...11. As mentioned in §3.3, the test is 
limited in the range p < 6. The error bars of Figure 7 are given by 
the maximum and minimum of the 80 subsamples. 

One can conclude that the fields of neutral hydrogen mass den- 
sity can be well described by the She-Leveque's scaling formula 
with statistical order as high as p = 6. It should be noted again. 
Figure 7 covers the physical scales from ~ 0.1 to 1.5 h~^Mpc. 
It indicates that cosmic neutral hydrogen fluid is turbulent on the 
scales considered. 



5 CONFRONTATION OF OBSERVED DATA WITH 
SIMULATION SIMPLES 

5.1 Lya spectrum synthesize 

Assuming ionization equilibrium under the uniform UV back- 
ground, we generate the field of neutral hydrogen fraction 
/hi(x) = nHi(x)/nH(x) at each cell, where nH(x) and nHi(x) 
are, respectively, the number densities of hydrogen and neutral hy- 
drogen at X. We synthesize 160 samples of normalized Lya trans- 
mitted flux ¥{z) = exp[— r(z)] from 2 = 4.9 to 2 = 6.3 using 
the same methods as Zhang et al (1997) and Paschos & Norman 
(2005). The optical depth r(z) is given by 

/oo 
nHi{x)V[z -X- v{x), b{x)]dx, (20) 
-00 

in which cro is the effective cross-section of the resonant absorp- 
tion and H is the Hubble constant at corresponding redshifts of 
the samples. The Voigt function is V[z — x — ^(a;), 6(a;)] = 
l/(7r^^^6) exp{— [2 — X — v{x)]^ /b^{x)}, where b{x) gives the 
thermal broadening. 

Along an randomly selected lines of sight, we synthesize an 
absorption spectrum F from 2 = 4.9 to 6.3 by dividing the spec- 
trum into redshift intervals of Az = 0.1. As the corresponding 
physical length scale for this redshift interval is larger than our sim- 
ulation box size, we integrate eq. (20) over the simulation dump 
periodically. Each spectrum is resolved with the same resolution 
of observation. Gaussian noise is added with signal-to-noise ratio 
S/N = 10. 

We calculate the Gunn-Peterson optical depth and its disper- 
sion of the 160 synthesized samples. The results are shown in Fig- 
ure 8. The observed results (Fan et al. 2006) are also shown in the 
Figure. The simulation samples basically are consistent with the 
observations except the dispersion at redshift 2 = 6.2. The devi- 
ation at redshift z = 6.2 might be due to the available observed 
samples at the redshift range z ~ 6.2 is too few. 

We also do the statistics of the probability distribution function 
of the transmitted flux ¥{z) of observation and simulation from 
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Figure 5. ln[_Fp+i(r)/F3(r)] vs. \n[Fp{r) / F2{r)] for simulation samples of the mass density field of neutral hydrogen at redshifts 5.0, 5.4, 5.8, 6.2. The 
statistical order to be p = 1, 1.5, 2, and 2.5. 
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mass density field. The number 7 is given in 



z — 5.3 to 2; — 6.1. The results are presented in Figure 9. The 
PDF of simulation samples can also fit observations within 1-sigma 
range. Thus, all low order statistics of the simulated samples are 
consistent with observation. 

5.2 Redshift dependence of /3 

Using the method of §3.2, we calculate the parameters /3 of simu- 
lated samples of Lya spectra. To mimic the observation samples. 



which contains totally 19 lines, we divide 160 one-dimensional 
simulation samples into 8 subsamples, each of which has 20 1-D 
samples. The uniform UV background eq.(17) and eq.(18) are used 
to produce the simulation samples. Figure 10 presents the redshift 
dependence of parameter /3 for observed data and simulation sam- 
ples. The error bars of the simulation results are given by the max- 
imum and minimum of the 8 subsamples. 

Although the redshift-dependence of the UV background 
eqs.(17) and (18) are different at 2; > 5, Figure 10 shows that 
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P=0.3 7=0.7 0=0.1 



P=0.31 7=0.65 0=0.22 - 



P=0.32 7=0.60 0=0.29 - 



P=0.35 7=0.60 0=0.26 - 



Figure 7. Intermittent exponent ^(p) of simulation samples of the mass density field of neutral hydrogen. The error bars are the variance of ^(p) over 80 
sub-samples, each of which contains of 2 one-dimensional samples. 
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Figure 8. Redshift-dependence of mean optical depth (left) and its variance (right). It shows observed result (circle) and simulation samples with UV back- 
ground eq.(17) (sohd Une) and eq.(18) (dotted Une). 



the redshift-dependence of /? given by the two models are about 
the same. It is because the uniform field J21 would not change the 
non-Gaussianity significantly. The simulation results of parameter 
/3 are basically consistent with the observation, only at z = 5.4 
shows a small deviation. Therefore, the ACDM universe embedded 
with a uniform UV background history is able to explain the feature 
of /3 redshift-revolution. 



Figure 7 shows that the parameters /3 of neutral hydrogen den- 
sity field is in the range 0.30 — 0.35, which is much less than the 
results shown in Figure 10. It indicates that the non-Gaussianity 
of neutral hydrogen field is different from that of Lya transmitted 
fiux. This is because the non-Gaussianity of Lya transmitted fiux 
depends not only on the mass density field of neutral hydrogen, but 
also on the velocity field via the Voigt convolution eq.(20). 



5.3 Parameter 7 

Similar to §3.3, we use high order moment {St^^) to quantify Kp 
(eq.(15)) and then the parameter 7. The results are listed in Table 
2. The values of 7 for p — 2 and 3 in each z range are the same 
within their errors. Comparing Table 1 and 2, the parameter 7 given 
by observation data and simulation sample are also consistent with 
each other within their errors. 

The parameter 7 for observational samples of Lya transmitted 
flux at « = 2.5 is found to be 0.58 ± 0.20 (Lu et al 2009), which is 
higher than that Usted in Table 2. It might indicate that parameter 7 
is decreasing with redshift. However, we should keep in mind that 
the error bars of observational sample are large both at high and low 
redshift, the redshift evolution of parameter 7 is not certain yet. 

Figure 6 reveals that parameter 7 of neutral hydrogen is in 
the range 0.80 — 0.89, which is much larger than that shown in 
Tables 1 and 2. Parameter 7 measures the singular structures and 
higher 7 represents stronger singularity (Liu & Fang, 2008; Lu et al 
2009). Therefore, the field of neutral hydrogen mass density con- 
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Figure 9. Probability distribution functions of the transmitted flux at 2 = 5.3 — 5.5 (top left), 2 = 5.5 — 5.7 (top right),2 = 5.7 — 5.9 (bottom left), 
z = 5.9 — 6.1 (bottom right). Cross are observation data. The solid and dashed line in the center are results of samples generated with eq.(17) and eq.(18) 
respectively. The above and below set of Unes gives the 1 it errors. 
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Figure 10. Parameter /3 vs. redshift for real data of Lya transmitted flux (§3.2)(circle); and hydrodynamic simulation results samples with UV background 
eq.(17) (square), and eq.(18) (trigon). 



tains much more singular structures than that of Lya transmitted 
flux. 



5.4 Spatially non-uniform UV background 

Whether the UV background at redshift z ^ 5 — 6 is spatially 
uniform is an important problem of the history of reionization. 
With only the redshift evolution of the optical depth and its vari- 
ance at a ~ 5 — 6, it seems to be difficult to distinguish models 
with uniform background from inhomogeneous one. We will ex- 
plore whether the log-Poisson non-Gaussian feature is dependent 
on the inhomogeneity of the UV background. 

Many works have been done on the fluctuation of hydrogen- 



ionizing radiation background field. Analytical method and Monte 
Carlo simulations with randomly distributed sources were firstly 
used to estimate the intensity fluctuations and its effect on the 
Lya forest (Zuo 1992; Fardal & ShuU 1993; Croft et al. 1999). 
Large scale N-body simulations and hydrodynamical simulations 
then were applied to study the fluctuations of the UV background 
(Gnedin & Hamilton 2002; Meiksin & White 2003; Meiksin & 
White 2004;Croft 2004; Bolton et al. 2006). Investigations have 
been carried out recently on this problem at high redshift (Wyithe & 
Loeb 2005; Mesinger & Furlanetto, 2009; Furlanetto & Mesinger 
2009). The general features of an inhomogeneous UV backgroimd 
can be sketched as 

J2l(x,«) = J2l(«)[H-<5(x,2)], (21) 
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Table 3. Parameter /3 and UV background models 



z 


A 


B 


C 


D 


5.4 


0.48± 0.02 


0.57± 0.02 


0.56± 0.02 


0.55± 0.02 


5.8 


0.52± 0.02 


0.50± 0.02 


0.53± 0.02 


0.52± 0.02 



z 


P 






5.0 


2 


-0.29± 0.04 


0.27; 




3 


-1.71±0.06 


0.29; 


5.4 


2 


-0.29± 0.06 


0.29; 




3 


-0.73± 0.12 


o.3o; 


5.8 


2 


-0.23± 0.03 


o.2o; 




3 


-0.59± 0.07 


0.23; 


6.2 


2 


-0.15±0.11 


0.14; 




3 


-0.38± 0.24 


0.15; 



Table 2. Parameter - at redshift = 5.0 to 6.2 



+0.04 
-0.04 
+0.02 
-0.03 
+0.06 
-O.O,"") 
+ 0.05 
-0.04 
+0.03 
-0.03 
+0.03 
-0.03 
+0.10 
-0.10 
+0.09 
-0.10 



where J2i{z) is the UV background given by eqs.(17) or (18). 
Equation (21) means, the field of the UV background is fluctuated 
with respect to its mean J21 (z). 

Unfortunately, no commonly accepted J21 (x, z) or (5(x, z) 
are available. The fluctuation field 5(x, z) may or may not be corre- 
lated with the density distribution of cosmic baryon matter. In this 
context, we consider 4 toy models of 5(x) as follows: 

• A. 5(x) = 0.1(;(x), where g(x) is a Gaussian random field 
with variance = 1.0. 

• B. (5(x) = 0.3(/(x). 

• C. (5(x) is given by 



5(x) = 



So sin(f min(lg(p(x))/2.0, 1.0)) p(x) > 1.0 
5o(f max(lg(p(x))/2.0, -1.0)) p(x) < 1.0 



(22) 



where So = 0.1. 
• D. The same as model C, but So = 0.3. 

For models A and B, the fluctuations of UV background are statis- 
tically independent of the IGM density field, while it is correlated 
to the IGM distribution for models C and D. 

Because the radiative transfer has not yet been included in our 
simulation, we simply add the non-uniform UV background to the 
outputs of simulations with eq.(17) and produce samples of the neu- 
tral hydrogen density field and Lya transmitted flux at redshift 5.4 
and 5.8. Since our goal is only to see whether the non-Gaussianity 
is affected by a fluctuating UV background, the post processing 
method would be acceptable as a first try. 

The j3 values of the Lya transmitted flux samples based on 
non-uniform UV background eq.(22) are listed in Table 3. Al- 
though the l3 values among models A,B,C and D sometimes show 
l-cT deviation from each other, these simplified models seem to 
have no effect in general. It may indicate that the effect of inhomo- 
geneous UV background on the log-Poisson non-Gaussian feature 
is not detectable yet with current data. 



trace the collisionless dark matter. The non-Gaussianity of the mass 
and velocity fields of baryon fluid is given by the hydrodynamics of 
the cosmic flow. A common property of the evolution of a Navier- 
Stokes fluid is to reach a turbulent state when the Reynolds num- 
ber is high (e.g. Zhu et al. 2010). In the scale free range, the fully 
developed turbulence is of statistically quasi-steady characterized 
by log-Poisson hierarchy. In this regime, the cosmic baryon fluid 
undergoes the evolution of clustering and finally falls into massive 
halos of dark matter to form structures, including light-emitting ob- 
jects. 

The cosmic baryon fluid or the IGM, both at high and low red- 
shifts, exhibits the log-Poisson non-Gaussian features in the range 
from the onset scale of the nonlinear evolution to the dissipation 
scale, i.e. the Jeans length. The log-Poisson non-Gaussianity has 
been identified in various simulated and observed samples of the 
IGM, including the mass density and velocity fields of baryonic 
matter, the density field of neutral hydrogen, and Lya transmitted 
flux. Although these fields are different from each other, they show 
the common behavior of log-Poisson non-Gaussianity. This sce- 
nario is further supported with high redshift Lya transmitted flux. 

The log-Poisson non-Gaussian parameter (5 of the mass den- 
sity field on the same physical scales is found to be increasing with 
redshift (Liu & Fang 2008). It implies that the mass density field 
of the IGM is less non-Gaussian at high redshifts. The log-Poisson 
non-Gaussianity at high redshifts should be weak than that of low 
redshifts on the same physical scale. However, we found in this pa- 
per that observed data of Lya transmitted flux at redshift 2 = 5-6 
shows about the same level of log-Poisson non-Gaussianity as low 
redshift z ^ 2 data. This is because the scales covered by the data 
at z = 5 - 6 are smaller than the data at ^; ~ 2. The weak evolu- 
tion of p given by this paper implies that turbulence state of cosmic 
baryon fluid at 2 ~ 2 on physical scales 1-10 h~^ Mpc is about 
the same as that at 2 = 5 - 6 on physical scales 0.1 - 1 h~^ Mpc. 
This is reasonable considering that the Jeans length and the typical 
scale of onset of nonlinear evolution are increasing with time. The 
log-Poisson non-Gaussianity at low and high redshifts are about the 
same once the nonlinear evolution is fully developed. 

In a word, on the scales of fully developed turbulence, the pa- 
rameter f5 of the IGM at high and low redshifts should be about the 
same. This property is very different from the Gunn-Peterson opti- 
cal depth and its variance, both show a strong redshift dependence. 
The Gunn-Peterson optical depth and its variance are sensitive to 
the intensity of the UV background. Oppositely, the parameters of 
the log-Poisson non-Gaussianity are weakly dependent on the in- 
tensity of the UV background. A more self-consistent and delicate 
handle of the inhomogeneity of the UV background in the cosmic 
hydrodynamic simulation and much more high quality observation 
data would help us to use the log-Poisson non-Gaussianity to in- 
vestigate the ionizing background at high redshift. 



6 DISCUSSION AND CONCLUSION 

Nonlinear evolution of baryon fluid at high redshift is a central 
problem of cosmology. It is well known that in the nonUnear regime 
the dynamical behavior of cosmic baryon fluid does not always 
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